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I.  INTRODUCTION 


Determination  of  the  physical  properties  of  projectiles  (mass,  center  of 
gravity  and  mass  moments  of  inertia)  is  required  in  each  of  the  particular 
fields  of  ballistics  design  *  Whether  the  projectile  is  in  the  bore  of  a  gun, 
in  flight,  or  at  the  target,  the  mechanical  performance  is  constrained  by  the 
limitations  imposed  by  the  physical  properties  •  Within  the  Ballistic  Research 
Laboratory  (BRL)  a  FORTRAN-based  calculation  of  these  critical  quantities  has 
been  employed  for  many  years,  but  the  advent  of  the  microcomputer,  with 
appropriate  programming,  now  provides  instant  access  to  computational 
facilities  in  BASIC*  This  report  presents  a  program  for  use  in  the  HP9845 
system  which  will  generate  values  for  the  physical  properties  of  typical 
projectiles . 


II.  PROCEDURE 

Conversion  of  the  Neitzel^  program  to  BASIC  consists  of  writing  BASIC 
statements  comparable  to  the  FORTRAN  equivalent  .  The  FORTRAN  "FUNCTION” 
statement  is  eliminated  through  the  use  of  subroutines  .  "DO"  loops  are 
replaced  by  "FOR* . . .NEXT”  loops.  "PRINT"  statements  are  substituted  for 
"WRITE"  commands .  The  scheme  of  dividing  the  projectile  into  sections  was 
changed.  The  projectile  is  segmented  along  the  horizontal  axis  to  provide 
boundaries  for  integration  and  the  Neitzel  program  used  the  initial  points 
from  each  input  region  as  boundary  points.  The  BASIC  version  differs  in  that 
the  division  into  segments  proceeds  from  the  lowest  to  the  highest,  these 
points  being  either  terminal  or  initial  points  and  creating  more  regions  . 
After  division  into  sections,  each  section  is  layered  and  integrated.  The 
methods  employed  are  identical  for  both  systems. 

Fins  are  approximated  as  equivalent  thin  plates  spaced  uniformly  around 
the  projectile  axis  of  rotation.  With  the  exception  of  the  transverse  mass 
moment,  any  number  of  fins  may  be  stipulated.  The  calculation  of  the 
transverse  mass  moment  of  the  fins  is  exact  for  either  4  or  6  fin  blades  . 

Data  is  entered  into  the  program  in  response  to  display  commands  .  Each 
segment  of  the  axisymmetric  projectile  is  identified  as  conical  (to  include 
cylinders  of  zero  slope)  or  toroidal  (circular-arc  cross  section),  fixed  by 
initial  and  terminal  location  and  assigned  a  material  density.  The  fin 
section  entry  considers  only  plates  of  flat  geometry.  A  fin  thickness  entry 
is  also  required.  Any  consistent  system  of  units  is  valid.  Before 
calculation  begins,  a  user  check  is  imposed.  All  input  data  is  displayed  and 
corrections  entered  as  required  without  rerunning  the  complete  program. 

III.  RESULTS 

The  results  are  illustrated  by  means  of  two  sample  runs.  Figure  1  is  a 
schematic  of  a  typical  projectile.  Figure  2  shows  an  outline  drawing  of  a 
typical  HEAT  round  and  Table  1-A  provides  a  complete  list  of  the  input  and 
output  data  using  this  program.  The  English  units  are  used  for  this  example 


^G.P.  Neitzel,  "A  Computer  Program  to  Calculate  the  Physical  Properties  of  a 
System  of  Co-Axial  Bodies  of  Revolution ,n  Ballistic  Research  Laboratories 
Memorandum  Report  No.  2215 >  August  1972  (AD  904Z78L) . 

5 


»  V 

r  " 


Figure  1.  Projectile  Calculating  Nomenclature 
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Figure  2.  120  mm  HEAT  Round 


TABLE  1-A.  HP9845  LIST  OF  INPUT  AND  OUTPUT  FOR  120  mm  HEAT  ROUND 
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to  allow  direct  comparison  with  the  furnished  reference  drawings.  Table  1-B  is 
the  output  for  the  same  data  using  the  program  from  Reference  1  and  the  CDC 
main  frame  computer  facility.  Figure  3  and  Table  2  illustrate  the  application 
to  a  representative  kinetic  energy  projectile  for  which  measured  values  are 
available.  The  results  are  given  in  both  metric  and  a  caliber^  notation. 
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W.F.  Donovan ,  " One  Factor  Affecting  the  Dispersion  of  Long  Rod  Penetrators ,  " 
ARBRL-MR-03020,  May  1980  (AD  086095) 
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Figure  3.  Line  Drawing  of  Flechette  Bullet  Projectile 
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APPENDIX  A 
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460 

470 

430 

490 

500 

510 

520 

530 

540 

550 

560 

570 

530 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

749 

750 
760 
770 


X2 


X2 


print  page 

PRINT  "  A0  A1  ** 

PRINT  USING  B*1 o| I , A0< I > » A1 < I > , X0< I > 

PRINT 

PRINT  "ENTER  TERMINAL  POINT" 

INPUT  Xf<I> 

PRINT  PACE 

PRINT  "  R0  A1  *1 

PRINT  USING  B*  1  o  J  I,A0(I)|Al(I)|X0<O|Xt  <I) 

PRINT 

PRINT  "ENTER  DENSITY" 

INPUT  Rod) 

PRINT  PAGE 

bdtWT1*  A0  ft  1  XI  *  *  * 

PRINT  USING  Btlo;  I,ft0<I),ftl<I>»X0<I:>>Xf<I)»Ro<I 
PRINT 

PRINT  " ENTER  FINAL* 

INPUT  Final 
PRINT  PAGE 
REM 

PRINT  "AND  WHETHER  THE  ARC  IS  CONCAVE  UP  <  INPUT  A  2>  OR 
PRINT  "CONCAVE  DOWN  (INPUT  A  1>.“ 

REM 


Ro" 


Ro" 


Ro" 


PRINT 

PRINT 


PRINT  “ENTER  X  COORDINATE  OF  THE  ARC  CENTER" 
INPUT  Xc  < I > 

PRINT  PAGE  y1 

PRINT  “  Xc  Yc  Rc  '  1 


it  y" 
780 
790 
300 
310 
820 
830 
i  ty" 
840 
850 
860 
870 
380 
890 
i  ty" 


PRINT  USING  Btlo; I,Xc<I> 
PRINT 

PRINT  "ENTER  Y  COORDINATE  OF 
INPUT  Yc<I> 

PRINT  PAGE 

PRINT  "  Xc  Ye 


THE  ARC  CENTER" 


Re  XI 


PRINT  USING  Bt 1 oj I , Xc < I > , Ye < I > 
PRINT 

PRINT  "ENTER  RADIUS  OF  ARC" 
INPUT  Rc(I) 

PRINT  PAGE 

PRINT  "  Xc  Ye 


XI 


Ro 


X2  Ro 


Cone iv 


Cone iv 


Cone  av 
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900  PRINT  USING  Be 1 oj I , Xc < I ) , Yc <  I  > , Rc < f>  " 

910  PRINT 

920  PRINT  "ENTER  INITIAL  POINT" 

930  INPUT  X0<I> 

940  PRINT  PAGE 

930  PRINT  ■  Xc  Yc  Rc  XI  X2  Ro 

i  ty" 

960  PRINT  USING  Btl o; I , Xc < I >, Yc < I > ,Rc < I > , X8< I > 

970  PRINT 

980  PRINT  "ENTER  TERMINAL  POINT" 

990  INPUT  Xf<I) 

1000  PRINT  PAGE 

1010  PRINT  "  Xc  Yc  Rc  *  XI  X2  Ro 

i  ty" 

1020  PRINT  USING  Be  1 o; I , Xc < I > , Yc < I> , Rc < I > , X0< I > , Xf < I > 

1030  PRINT 

1040  PRINT  "ENTER  DENSITY" 

1050  INPUT  Ro < I > 

1060  PRINT  PAGE 

1070  PRINT  "  Xc  Yc  Rc  XI  X2  Ro 

i  t  y  " 

1080  PRINT  USING  Bt 1 o; I f Xc < I > , Yc < I > , Rc < I > , X0< I > , Xf < I > , Ro< I > 

1090  PRINT 

1100  PRINT  "ENTER  CONCAVITY" 

1110  INPUT  Cone av<  I ) 

1120  PRINT  PAGE 

1130  PRINT  "  Xc  Yc  Rc  XI  X2  Ro 

i  ty" 

1140  PRINT  USING  Bel o; I , Xc <1 > , Yc < I > , Rc < I > , X0< I > , Xf ( I > , Ro< I ) , ConcavC I > 

1130  PRINT 

1160  PRINT  “ENTER  FINAL" 

1170  INPUT  Final 

1180  PRINT  PAGE 

1190  A0U)«0 

1200  At  < I )«0 

1210  PRINT  PAGE 

1220  IF  Fin*l-1  THEN  1248 

1230  NEXT  I 

1231  IF  Final-0  THEN  200 

1240  PRINT  “A  LIST  OF  THE  INPUT  VALUES  WILL  APPEAR  ON  THE  SCREEN  WHEN  YOU 
1250  PRINT  "PRESS  CONTINUE.  AFTER  REVIEWING  THE  DATA,  INDICATE  IF  YOU  WISH 
1260  PRINT  "TO  MAKE  ANY  CHANGES. <Y  OR  N>" 

1278  PAUSE 
1280  PRINT 

1290  PRINT  "  A0  Ai  XC  YC  R  XI  X2 

DENSITY" 


*  * 


Cone  ao 


Cone  a Kt 


Cone  au 


Cone  ao 
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V 


,  ajuu.  coil  ,  sx,  odd.  odd 

1330  NEXT  G 


2X, DDD . DDD  , 2X, DDD . DDD 
0  i 4X} DDD. DDD, 


1331  PRINT 

1332  PRINT 

1333  PRINT  "PLEASE  INDICATE  IF  YOU  WISH  TO  MAKE  ANY  CHANGES.  <Y  OR  N)" 
1340  INPUT  Chang## 

1330  IF  Changt#* "N"  THEN  1510 

1331  PRINT  "ENTER  INDEX,  NAME  AND  NEW  VALUE  OF  PARAMETER  TO  BE  CHANGED" 

1352  PRINT 

1353  INPUT  In, Na»t$, Val u# 

1354  IF  Namt#<  > " A0"  THEN  135? 

1 335  A0< I n ) *Val ut 

1336  GOTO  1290 

1357  IF  Namt#<  >"Al "  THEN  1360 
1338  A1 < In>*Va1 ut 

1359  GOTO  1298 

1360  IF  Namt$<  >"X1 "  THEN  1363 

1361  X0  < I n ) *Val ut 

1362  GOTO  1290 

1363  IF  Namt $<  > " X2"  THEN  1366 

1364  Xf < I n) *Val ut 

1365  GOTO  1290 

1366  IF  Namt$<  > "Ro "  THEN  1373 

1 367  Ro  < I n ) *Val ut 

1368  GOTO  1290 

1373  IF  Namt$<  >"XC"  THEN  1376 

1374  Xc  < In)«Va1 ut 

1375  GOTO  1290 

1376  IF  Namt$<  >" YC"  THEN  1379 

1377  Yc  < In )«Val ut 

1378  GOTO  1290 

1379  IF  Namt#<  > "R"  THEN  1382 

1380  Rc<In)*Valut 

1381  GOTO  1290 

1382  IF  Namt #<  > " X 1 "  THEN  1385 

1 383  X0  < In) »Val ut 

1384  GOTO  1290 

1385  IF  N am t # <  > " X 2 "  THEN  1388 

1386  XT < I n >*Va1 ut 

1387  IF  Namt#<  > "Ro H  THEN  1390 

1388  Ro(In)*Valut 

1389  GOTO  1290 

1390  IF  Namt$<  >"Conc au"  THEN  1290 
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1391  Conc*M< In)*Val ue 

1392  GOTO  1290 

1393  IF  ConcauCLi ne  >*0  THEM  1460 

1400  PRINT  “PLEASE  INPUT  THE  X- , Y-COORD INATES  OF  THE  CENTER  OF  THE  ARC," 

1410  PRINT  “THE  RADIUS  OF  THE  ARC,  INITIAL  PT.  , TERMINAL  PT. , THE  DENSITY  “ 

1420  PRINT  "AND  WHETHER  THE  ARC  IS  CONCAVE  UP  (INPUT  A  2>  OR" 

1430  PRINT  "CONCAVE  DOWN  (INPUT  A  1>." 

1440  INPUT  Xc (Li ne >,Yc(Line),Rc(Line),X0 (Line), Xf (Line), Ro (Line), Cone ao (Line) 
1450  GOTO  1480 

1460  PRINT  "PLEASE  INPUT  THE  Y- I NTERCEPT , SLOPE , I N I T I AL  PT., TERMINAL  PT . , &DENS I TY 

U 


1470  INPUT  A0(Line), Rl(Line>,X0(Line),Xf(Une),Ro(Line) 

1480  PRINT  "WOULD  YOU  LIKE  TO  MAKE  ANY  OTHER  CHANGES? (Y  OR  N>" 

1490  INPUT  Maybe* 

1500  IF  Maybe**" Y"  THEN  1370 
1510  Ti*I 
1520  PRINT 
1530  N® I 

1540  PRINT  "WOULD  YOU  LIKE  A  COPY  OF  THE  INPUT  DATA?( Y  OR  N)" 

1550  INPUT  Data* 

1560  PRINT  PAGE 

1570  IF  Dat a*® "N"  THEN  16O0 

1580  PRINTER  IS  0 

1590  PRINT 

1600  PRINT  "  A0  A1  XC  YC  R  XI  X2 

DENSITY" 

1610  PRINT 
1620  Li*N 
1630  FOR  X*1  TO  N 
1640  X(I)*X0(I) 

1650  J-N+I 
1660  X(J)«Xf(I) 

1670  NEXT  I 
1680  Nn*2*N 
1690  N*Nn 
1700  L-l 

1710  FOR  1*1  TO  N 
1720  F 1 op*0 
1730  GOSUB  P»in 
1740  GOTO  1930 
1750  Pmi  n : REM 
1760  Jstart*! 

1770  Xm i n® 1 . 0E99 
1780  FOR  V* J*t art  TO  Ti 
1790  IF  X( V ) >*Xmi n  THEN  1830 
1800  Xm  i  n*Xj^V ) _ 
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1810  Indt x«V 

1820  Q-V+Tt 

1830  NEXT  V 

1840  X*ave*X< Jstart > 

1830  X2*ave«X< Jetart+TI ) 

1860  X3s4u*"A0< J*t art ) 

1870  X4*ave*ft 1 < Jst art > 

1880  X3*ave«Xc < Jet 4rt > 

1890  X6*aoe»Yc < Jetart ) 

1900  X7*ave *Rc < Jst art > 

1918  X8»aye*Ro< Jst art ) 

1920  X9*ave*Concao< Jet 4rt ) 

1930  X<J*tart)*X<Index) 

1948  X< Jstart+Ti )*X<Q) 

1930  fl0< Jst 4rt ) »fl0< Index) 

1960  ftl < Jet  art )«ftl < Index) 

1970  Xc <Jst4ri )«Xc < Indtx) 

1980  Yc <J*tart )*Yc < Index) 

1990  Re < Jet  art )*Rc < Index) 

2000  RoCJstart )*Ro< Index) 

2010  Cone au (Jet  art  )*Coneav<  Index) 

2020  X<  Index)*Xsa'y»e 
2038  X<Q)*X2eave 
2040  ft0< Index)*X3save 
2030  ft  1 < I ndex )®X4save 
2060  Xc < I ndex)»X3save 
2070  Yc < Index)*X6e ave 
2088  Rc(Index)*X7*aue 
2098  Ro< Index)*X8sa^e 
2100  Cone ao< Index >*X9eaoe 
2110  Xbp<L)«X< Jstart ) 

2120  Xfp<L)sX(Jstart+Ti > 

2130  Jet art*Jet art +1 
2140  L*L+ 1 

2150  IF  Jet  art  <  *T i  THEN  1770 
2160  FOR  G* 1  TO  Ti 

2170  PRINT  USING  Beli w; R0<G> , ft! <G) , Xc <G> , Yc <G) , Rc <C> , X<G> , Xf p<G> , Ro(G) 

2188  Bel i u : I MflGE  DDD. DDDD  , 2X, DDD. DDDD  , 2X, DDD. DDDD  , 2X , DDD. DDDD  , 2X , DBD . DDDD 

,  2X , DBD. DDDD  ,2X, DDD. DDDD  ,4X, DDD. DDDD 

2190  NEXT  G 

2200  PRINTER  IS  16 

2210  PRINT 

2220  PRINT 

2230  PRINT 

2240  PRINT  "THE  CftLCULftTIONS  ftRE  BEING  PERFORMED.  PLEftSE  EE  PftTIENT" 

2250  PRINT 

2260  Zstart*! _ _ _  _  _  _ _ _ _ 


*  v 


/ 
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2270  Xm i n2* 1 « 0E99 

2280  FOR  Z*Zstart  TO  Ti 

2290  IF  Xf p<Z> >*Xm i n2  THEN  2320 

2300  Xm i n2*Xf p  <  Z ) 

2310  Zoro*Z 
2320  NEXT  Z 
2330  Zs*ue*Xfp<Zst*rt > 

2340  XfpCZst4.rt)«Xfp<2oro> 

2330  Xfp<Zoro)*Zs*oe 

2360  Zst art *2st  *rt  +  1 

2370  IF  2st art <*T i  THEN  2270 

2380  Z*1 

2390  Count*0 

2400  FOR  Order* 1  TO  Ti 

2410  IF  Xfp<Z)-Xfp<Z+l>  THEN  2488 

2420  IF  Order* 1  THEN  2430 

2430  IF  X  <  Order  )*Cm»i<  Count  )  THEN  2310 

2440  IF  X<Order>*Xfp<Z>  THEN  2510 

2430  Count “Count ♦ 1 

2468  IF  X<  Order  ><Xfp<Z>  THEN  2300 

2470  Cmm< Count >«Xfp<Z> 

2480  Z«Z+1 

2490  GOTO  2410 

2508  Cmm< Count ) »X<Order ) 

2510  NEXT  Order 

2520  IF  Xfp<Z>>*X<Ti >  THEN  2570 
2330  IF  Z*Ti+l  THEN  2610 
2540  CmmCCount >*Xf p<Z> 

2350  IF  Xfp<Z)*Xfp<Z-l >  THEN  2390 

2360  Z*Z+1 

2370  Count *Count ♦ 1 

2580  GOTO  2330 

2390  2*Z+ 1 

2600  GOTO  2530 

2610  L*T i ♦ 1 

2620  X0<l>*Xbp<n 

2630  Xf <t)*Xfp<L> 

2640  1*1 
2630  I i *  1 

2660  Xt<  I  >*<Cmm<I)+Cmm<I  +  l>>/'2 

2670  K* 1 

2680  Cue* l 

2690  FOR  J«l  TO  Li 

2700  IF  <X<J>>Cmm< I))  OR  <Cmm< 1  +  1 > >X< J+Ti > >  THEN  2800 
2710  IF  Rc  < J)*0  THEN  2770 
2720  IF  Cone *u<  J  >»2  THEN  2730 
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2730  Fy<K>iYc<J>VSfilR<Rc<J>A2-<Xt < I >-Xc < J> )A2)+ft0< < J)*Xi < I > 
2740  GOTO  2780 

2730  FyOO«Yc<J)-SQR<Rc(J>/'k2-<Xt<I)~Xc<J>)A2>  +  fi0<J)+ftl(J)*Xt<I> 

2760  GOTO  2780 

2770  Fy <K)»A0< J)+fll <J)*Xt  <  I  > 

2780  M n  < K ) ■ J 

2790  K*K+1 

28O0  NEXT  J 

2810  K«K-1 

2820  JJ-1 

2830  J»1 

2840  X«*x*-2E99 

2830  FOR  R-J  TO  K 

2860  IF  Fy <R><Xm*x  THEN  2890 

2870  Xm*x-Fy<R> 

2880  RUh»R 
2890  NEXT  R 
2900  R**vt-Fy<J> 

2910  Fy<J)»Fy<Ri ch) 

2920  Fy<R1ch>*R 
2930 

2940  Mn< J>*Mn<Ri ch> 

2950  Mn<R i ch 
2960  J«J+1 

2970  IF  J<«K  THEN  2840 
2980  Guo*0 
2990  B«t  4*K 

3000  IF  B*t*«l  THEN  3050 
3010  Quo*Quo+l 
3020  Nup*Mn<Quo> 

3030  N 1 o»Mn ( Quo* 1 ) 

3040  GOTO  3070 
3030  Nup*Mn<Quo+l ) 

3060  N1 o*0 
3070  XI <I>«C»m<I) 

3080  Xu<I>«Cm»<I+i; 

3090  FI i p*0 
3100  GOSUB  Rfflbg \ n 
3110  FI ip«l 

3120  Xm»Xm  +  Pt  *Ro<Nup)*Fi 

3130  Bon*Fi 

3140  GOSUB  Rwbgin 

3150  F ] i p*2 

3160  Xcg*Fi/Bon 

3170  GOSUB  Rmbgin 

3180  FI ip«3 
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3190  Ami »A« i ♦. 5*Pi *Ro<Nup)*Fi 

3200  G i * F i 

3210  GOSUB  Rmbgin 

3220  2f*2f +Bon*Xcg*Pi *Ro<Nup> 

3230  Bb*Bb+Pi *Ro<Nup)*< . 23*Gi +F i > 

3240  B*ta*Btta-l 

3230  IF  Xu <I)*Xfp< Count-1)  THEN  3280 

3268  IF  Beta»0  THEN  3300 

3270  GOTO  3290 

3280  IF  N 1 o*0  THEN  3330 

3290  GOTO  3000 

3300  I«I+l 

3310  IF  I *L  THEN  3330 
3320  GOTO  2650 
3338  Cgproj*2f/Xm 
3340  B*Bb-Xm*Cgpro j A2 
3330  PRINT  PAGE 
3360  PRINT 

3370  PRINT  “DOES  THIS  BODY  HAVE  FINS  WHICH  NEED  TO  BE  INCLUDED  “ 

3380  PRINT  "IN  THE  CALCULAT I ONS?< Y  OR  N>" 

3390  INPUT  Contn* 

3391  IF  Cont n$*M Y“  THEN  Data*-"Y" 

3400  PRINT  PAGE 

3410  REM 
3420  PRINT 

3430  PRINT  H  DO  YOU  WANT  A  HARD  COPY  OF  THE  RESULTS  <Y  OR  N>" 

3440  INPUT  Out $ 

3450  IF  Out  $*" Y"  THEN  PRINTER  IS  0 

3460  PRINT  -  MASS  CG  AX  MOM  TRANS  MOM 

3470  PRINT  USING  Format ; Xm , Cgproj , Rm i , B 

3480  Format : IMAGE  6X , DDDD . DDDD , 3X , DDDD. DDDD, 5X, DDDD. DDDD , 5X , DDDD. DDDD 

3481  IF  Contn#»"Y"  THEN  4360 

3482  IF  Cont n$*HNM  THEN  6190 
3490  PRINTER  IS  16 

3300  PRINT 
3510  PRINT 

3520  PRINT  "IF  YOU  WISH  TO  START  AGAIN,  PRESS  RUN." 

3530  END 

3340  Rmbgin:  DIM  AC9),B<9) 

3530  FOR  M-l  TO  9 
3360  A<M) *0 
3370  B < M) »0 
3380  NEXT  M 
3399  REM 
3600  X1*X1 (I) 

3610  X*X1  _ _ _  .  _ 
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3620  GOSUB  Fxx 

3630  GOTO  3890 

3640  FxxlREM 

3630  Xxc < Nup >*X-Xc <Nup) 

3660  IF  Rc  <Nup><  >0  THEN  3680 

3670  Xxc<Nup)*0 

3680  Xxc<Nlo)*X-Xc<N1o) 

3690  IF  Rc (N1 o)<  >0  THEN  3710 
3700  Xxc  <N I o > *0 

3710  IF  Cone *u<Nup ) *2  THEN  3760 

3720  Yu*(Yc<Nup)+<Rc<Nup)A2-Xxc  <Nup>A2) A.  5+R0<Nup>+ftl(Nup)*X)A2 
3730  IF  Cone *u<N 1 o  >  *2  THEN  3780 

3740  Y1 *<Yc<N1 o>  +  <Rc  <N1 o> A2-Xxc  <N1 o>a2> a. 5+R8<N} o>+Rl <N] o>#X> a2 
3730  GOTO  3790 

3760  Yu*  <  Ye  <  Nup  )-(Rc<Nup)A2-Xxc  (Nup) A2>A.  3+R0<Nup>+fil  <Nup>*X>A2 
3770  GOTO  3730 

3780  Y1 «<Yc<N1 o>-<Rc  <N1 o)A2-Xxc<N1o>A2)A. 3+R0<NJ  o>+fll <N1 o>*X>A2 

3790  Fx*Yu-Y 1 

3800  IF  FI ip*0  THEN  3880 

3810  IF  FI  ip-1  THEN  3870 

3820  IF  FI ip«2  THEN  3850 

3830  Fx-Fx*XA2 

3840  GOTO  3880 

3830  Fx*Yua2-Y1a2 

3868  GOTO  3880 

3870  Fx*X#Fx 

3880  RETURN 

3898  Fa-Fx 

3980  U1 -Xu  <  I ) 

3910  X*U 1 

3920  GOSUB  Fxx 

3930  F«Fx 

3940  H-U1-X1 

3930  R  < 1 ) * . 3*H*<Fa+F  > 

3960  Ip*l 

3970  Ic *8 

3980  !s*l 

3998  Ic - 1 

4000  H 1 *H 

4010  H*. 5*H 

4020  X-X1+H 

4030  Sum-0 

4040  FOR  Q* l  TO  I* 

4030  GOSUB  Fxx 
4060  $um*Fx+Sum 
4070  X-Hl+X 
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4088  NEXT  Q 
4890  I*=Is+I* 

4100  B<l>-.5*<A(l)+Hl*Su«i> 

4110  C*4 

4120  FOR  J-l  TO  Ip 
4130  K* J+  1 

4140  BOO-<C*B<J>-A<J>)/<C-1) 

4130  C«4*C 
4160  NEXT  J 
4170  FOR  J-l  TO  Ip 
4180  K-J+l 

4190  Abc*ABS<  <B<J)-B<K>>/BOO> 

4200  To  1  - 1 . 0E-5 

4210  IF  Abe -Tol <0  THEN  4330 

4220  Abc*AB$<<AOO-B<IO>/B<K>> 

4230  IF  Abc-Tol <0  THEN  4330 
4240  NEXT  J 

4230  IF  Ip-8  THEN  4270 
4260  Ip*Ip+l 
4270  Ic-Ic+1 
4280  FOR  J«1  TO  9 
4290  A< J)«B< J> 

4300  NEXT  J 

4310  IF  Ic<10  THEN  4000 

4320  PRINT  “Rmbgin  DID  NOT  CONVERGE  IN  TEN  STEPS." 

4330  F i *B<  K ) 

4340  REM 
4338  RETURN 
4360  PRINTER  IS  0 
4370  PRINT 
4380  PRINT 

4390  IF  Cont " Y "  THEN  Dat*f*MY" 

4400  PRINT  “THE  RESULTS  FOR  THE  AXISYMMETRIC  PROJECTILE  ALONE  ARE  GIVEN  BELOW," 

4410  PRINT  “THE  FINS  ARE  CALCULATED  SEPARATELY  AND  THEN  ADDED  TO  THE  BODY" 

4411  PRINT 

4412  PRINT  “  A0  A1  XC  YC  R  XI  X2  R 

o  CONC “ 

4413  FOR  G- 1  TO  1-1 

4414  PRINT  USING  Be  1 u; G, A0<G> , Al <G> , Xc <G> , Ye <G> , Re <G> , X0(G> , Xf p<G) , Ro<G>  ,  Conciv< 
G) 

4415  Be  1 u :  IMAGE  DD, DDD. DDD, 2X, DDD. DDD, 2X, DDD. DDD, 2X, DDD. DDD, 2X, DDD. DDD, 2X, DDD. 
DDD,2X,DDD.DDD,2X,DDD.DDD,5X,D 

4416  NEXT  C 

4417  PRINT 

4420  PRINT  “PROJTILE  MASS  CG  AXMOM  TRANSMOM'* 

4430  PRINT  USING  Form# t ; Xm , Cgpro j t Am i , B  _ _  _  
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5350 

PRINT  “ENTER  INDEX, 

NAME 

5360 

PRINT 

5370 

INPUT  Zn,N4M«# 

,V*1ut 

5380 

IF  N*w#*O,,A0,‘ 

THEN 

5410 

5390 

A0<In>-V*1uf 

5400 

GOTO  5230 

5410 

IF  NtftttO'Al" 

THEN 

3440 

5420 

Al<In>-V*1ut 

5430 

GOTO  5250 

5440 

IF  N*m#l<  > “X 1 “ 

THEN 

3470 

5430 

X0< I n )*V4l ut 

5460 

GOTO  3230 

5470 

IF  N4*t*<>“X2“ 

THEN 

5300 

5480 

Xf CIn>»V4lut 

5490 

GOTO  5250 

5500 

IF  N*n«$<  > “Th“ 

THEN 

3530 

5310 

Th< In)»V4l ut 

5320 

GOTO  5230 

5530 

IF  N4mt«<>“Ro“ 

THEN 

5230 

3340 

Ro<In)«V4lut 

3350 

GOTO  5250 

5360 

PRINT 

3370 

PRINTER  IS  0 

3380 

M 4*4*0 

3590 

REM 

5600 

PRINT  "  A0 

Ai 

5610 

FOR  G* 1  TO  J 

5620 

REM 

3630 

REM 

3640 

REM 

3630 

PRINT  USING  Bt 1 o; G, A0<G> 

5660 

Bt 1 o; IMAGE  DD,: 

DDD. DDD 

XI 


X2 


Th 


Ro" 


IDD.DDD  ,2X,DDD.DDD  , 2X, DDD . DDD  , 2X, DDD. DDD 

t 2X, DDD . DDD  ,4X,DDD.DDD 

5670  IMAGE  2X, DDD. DDD  , 2X, DDD . DDD  , 2X, DDD. DDD  ,2X,DDD.DDD  , 4X, DDD. DDD 

5671  REM 
5680  REh 

5690  M4***Ro<G>*Th<G>*<Al <G>/2*<Xf <G)A2-X0<G> A2>+A0<G>*<Xf <G>-X0CG>  >  >+M4** 

5700  M4**3*M4ss*L 

5710  P-Ro<G>*Th<G>*<Al <G)/3*<Xf <G>A3-X8<G>A3>+A0<G>/2*<Xf (G>A2-X0<G)A2>  >  +  P 
5720  Cent gr*v»P/N*ss. 

5730  REM 
5740  REM 
5750  REM 
5760  REM 

5770  Ql*Ro<G)*Th<G> 

J5788  Q2*Al<G)A3<''12*<Xf<G>A4-X0<G)A4)+fll<G>A2*A0<G)/3*<Xf<G>A3-X0(G>A3) 
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U3*A1  <G)*<fi0<G>^2/'2+Th<G)-A>2/24>*<Xf  <G) A2-X0<G>A2> 

3800  Q4-R0<G)*<A0<:G>A2/3+Th<G:>A2/'12>*<Xf  <G>-X0<G>  > 

5810  Q5«Q1*<Q2+Q3+Q4>+Q5 
3820  Q*L*Q5 
3830  Axmom*Q 

'  3840  IF  L*4  THEN  3870 

5830  IF  L»6  THEN  5912 

5851  PRINT  "IF  L  IS  OTHER  THAN  4  OR  6,  THE  CALCULATIONS  ARE  NOT  EXACT  FOR  THE  TR 
\  ANSVERSE  MOMENT.  THE  CALCULATIONS  FOR  THE  OTHER  PROPERTIES  ARE  EXACT" 

^  5852  REM 

5860  REM 

5870  Ul«2*Ro<G>#Th<G> 

3880  W2»<A1 CG)^3^12+A1 <G>/2)*<Xf <G>a4-X0<G>a4>-K2/3*A0<C>+A1 <G > A2*A0< G>/3 ) *<Xf <G 
>A3-X0<G)A3> 

5890  U3-CA1  <G>#Th<G)  A2>'24+Al<G>*A0<G>A2/2>*<Xf  <G>A2-X0<G>A2> 

5900  W4* < A0 CG>^Th< G) ^2^1 2+A0CG) A3/3>*<Xf <t)-X0<G ) > 

5910  W3«W1*<W2+W3+W4>+W3 

5911  IF  L-4  THEN  3920 
3912  Wl»2*Ro<G)*Th<G) 

5913  W2*A1 <G)/4*3*<Xf (G) A4-X0<G) A4)+<A0<G>+fil <G>*Th<G>*<Xf <G>A3-X0<G) A3> ) 

5914  U3-<Al<G>/8*Th<G>A2+3*Al <G>*A0<G>*Th<G> >*<Xf <G)A2-X0<G> A2> 

5913  W4»<A0<G>*Th<G)A2/4+3*A0<G>A2*Th<G> >*<Xf <G)-X8<C>  > 

5916  W5»W 1 *< W2+W3+W4 >+W5 
3920  REM 
5930  REM 

5940  W6»W5*L-M*ss5*<P/M*ss)A2 

5930  Tr*nsmom*U6 

5960  REM 

3970  FIXED  4 

3980  NEXT  G 

5990  REM 

6000  PRINT 

6010  PRINT  "FINS  MASS  CG  AXMOM  TRANSMOM 

NUMBER" 

6020  PRINT  USING  Formi  t ;  M4.ss3,  Centgr-av,  Axmom,  Tra.nsmom,  L 

6030  For mi t: IMAGE  X, 5X, DDDD. DDDD, 5X, DDDD. DDDD, 5X, DDDD. DDDD, 5X, DDDD. DDDD, 14X.D 

6040  REM 

6050  REM 

6060  REM 

6070  REM 

6080  PRINTER  IS  16 
6090  REM 

6100  FOR  G»1  TO  J 
6110  REM 

6120  PRINT  USING  Be  1 o; G, A0<G) , A1 <G> , X0<G> , Xf <G> , Th<G> , Ro<G) 

6130  NEXT  G  _ 
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6131  REM 

6132  PRINT 

6140  PRINT  "Numbtr  of  fin*"(L 
6150  PRINT  “Mas*5" , M***5 
6160  PRINT  "Ctntgrav" , Ctnt gr*v 
6170  PRINT  "Axmom" , Axmom 
6180  PRINT  "Tr*nsmom“ , Trantmom 

6190  PRINTER  IS  0 

6191  PRINT 

6192  PRINT 

6208  Tot «]  ma*s«Ma**5+Xm 

6210  PRINT  "Tot *1  mass"  ,  Total  mass 

6220  Ntt  cg*(Xm*Cgpro  j+Mass5*Ctntgrav>/'Tot  al  m&ss 

6230  PRINT  "Nttcg", Nttcg 

6240  Tot al axmom»ftmi tAxmom 

6250  PRINT  " Tot 4l axmom M , Tot *1 axmom 

6260  Tot4Hr4n*mom«B+W6  +  M4**5*<Cgprog-C*ntgr4V>''2 

6270  PRINT  "Tot 4ltr*nj«o#" , Total transmom 

6271  PRINT 

6272  PRINT 

6280  PRINTER  IS  16 

"ENTER  THE  REFERENCE  DIAMETER  USING  THE  UNITS  EMPLOYED  IN  THE  TEXT" 
"  D i 4" , D 


"ENTER  THE  REFERENCE  SPECIFIC  WT  OF  WATER  USING  THE  UNITS  EMPLOYED  IN 
"Gamma" , G 


6380  Mass7*Tot al  mas */DA3/G 

6390  Ctnt gra^7=Ntt cg/D 

6480  Ax(aom7»Tot4l  4xmom/D/s5/G 

6410  Tr4n»mom7*Tot4ltr4nsmof»/D/s5/'G 

6420  PRINT 

6430  PRINTER  IS  0 

6440  PRINT  "M4S*7M,M4J*7> "C4l^3" 

6450  PRINT  "Ctnt grav?" , Ctnt grav? ,  "c 4l  * 
6460  PRINT  " Axmom?"  ,  Ax»om7,  "c4l  A5" 

6470  PRINT  "Transmom?" ,  Trans  mom  7,  "cal  '"S" 

6471  PRINT  "Dia",D,"cal  rtf" 

6480  PRINTER  IS  16 

6490  END 


6281  PRINT 
6290  PRINT 
6300  PRINT 
6310  INPUT 
6320  PRINT 
6330  PRINT 
6340  PRINT 
THE  TEXT" 
6350  INPUT 
6360  PRINT 
6370  PRINT 
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LIST  OF  SYMBOLS 


A#  intercept  on  vertical  axis 
A1  slope  of  boundary  segment 
L  number  of  fins 
R  radius 
Ro  density 

Th  average  fin  thickness 

XC  center  of  radius  along  horizontal  axis 

XI  initial  segment  reference 

X2  terminal  segment  reference 

YC  center  of  radius  along  vertical  axis 


A 


% 

N. 
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